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Abstract — Ant Colony Optimization (ACO) is a very popular 
metaheuristic for solving computationally hard combinatorial 
optimization problems. Runtime analysis of ACO with respect 
to various pseudo-boolean functions and different graph based 
combinatorial optimization problems has been taken up in recent 
years. In this paper, we investigate the runtime behavior of an 
MMAS*(Max-Min Ant System) ACO algorithm on some well 
known hypergraph covering problems that are NP-Hard. In 
particular, we have addressed the Minimum Edge Cover problem, 
the Minimum Vertex Cover problem and the Maximum Weak- 
Independent Set problem. The influence of pheromone values 
and heuristic information on the running time is analysed. The 
results indicate that the heuristic information has greater impact 
towards improving the expected optimization time as compared 
to pheromone values. For certain instances of hypergraphs, we 
show that the MM AS* algorithm gives a constant order expected 
optimization time when the dominance of heuristic information 
is suitably increased. 

Index Terms — Ant Colony Optimization; Hypergraph Set 
Cover; Runtime analysis. 

I. Introduction 

Ant colony Optimization (ACO), first proposed by Dorigo 
[1], is a class of nature-inspired stochastic algorithm which 
is widely used for various combinatorial optimization prob- 
lems. The ACO algorithm is derived from the food searching 
behaviour of real ant colonies, who have demonstrated their 
ability to find the shortest path from their nest to the food 
source by stigmergic communication via pheromones. For a 
given problem, the ACO algorithm constructs a solution by 
random walk on a construction graph, with the walk influenced 
by pheromone values and heuristic information. Thus, ACO is 
able to incorporate problem specific knowledge in terms of 
heuristic information which separates it from other classes 
of Evolutionary Algorithms (EAs), etc. Although ACO is 
not evolutionary in nature, it has by and large attracted the 
attention of Evolutionary Computing researchers because it 
is also a random search heuristic and both have overlapping 
application domains. Many different variants of the original 
ACO algorithm have been proposed and applied successfully 
to many real-world problems [2]. 

Initial theoretical studies on ACO focussed on convergence 
properties and speed of different variants of ACO. Stiitzle 
and Dorigo presented the convergence proof for the Max-Min 



ant system (MM AS [5]) in [4], while Gutjahr showed it for 
the Graph based ant system (GABS) [6]. The foundations of 
runtime analysis of ACO were laid in [8] [9]. The analysis 
is typically done in a fashion that has been widely used for 
randomized algorithms and for simple evolutionary algorithms 
like (1h-1) EA. In this approach, either the expected optimiza- 
tion time and/or the success probability after a certain number 
of steps are analysed. 

Investigations on the runtime performance of ACO first 
started for various pseudo-boolean functions [3,10,11,12] 
where a simple ACO variant 1-ANT was analysed like (1h-1) 
EA [7]. Another ACO variant MM AS has also been widely 
investigated for which the phase transition present in 1-ANT 
due to the evaporation factor does not occur [13,14]. The above 
techniques have been extended to different polynomial time 
combinatorial optimization problems, namely the Minimum 
Spanning Tree problem[15]. Single Destination Shortest Path 
problem[16,17] and Minimum-cut problem[18]. Among NP- 
Hard problems, certain instances of the Travelling Sales- 
man Problem have been investigated in [19,20]. A detailed 
overview of runtime behaviour of Randomized Search Heuris- 
tics (RSHs), (1h-1) EA and ACO on different pseudo-boolean 
functions can be found in [21] and on several graph based 
combinatorial optimization problems can be found in [22]. 

In this paper, we undertake the runtime analysis of ACO 
for covering problems on hypergraphs. A hypergraph is a 
generalization of a graph where an edge, known as hyperedge 
consists of a non-empty set of vertices with any number of 
elements. On a hypergraph, different covering problems are 
the vertex cover and the edge cover problems, both of which 
can be weighted or unweighted. Independent set problems 
are also related to covering problems. The decision versions 
of these problems are NP-complete, while the optimization 
versions (minimum vertex/edge cover/maximum independent 
set) are NP-Hard[29]. The vertex cover problem is also known 
as transversal or hitting set problem, where the edge cover 
problem is a special case of set cover problem for hypergraphs. 
Computing the transversal hypergraph has many applications 
in machine learning, game theory, indexing of databases, SAT 
problem, data mining and optimization. 

Several interesting results have been obtained by applying 
various approximation algorithms for determining the vertex 



cover and set cover problems on a hypergraph or special 
cases of hypergraphs. Khot has published several results on 
the approximability and hardness of vertex cover problem on 
k-uniform hypergraphs[24,26]. In [24], the hypergraph vertex 
cover problem is shown to have apphed directly to several 
scheduUng problems : concurrent open shop and makespan 
minimization of assembly line. Okun[25] discusses the ap- 
proximation of hypergraph vertex cover for hypergraphs of 
bounded degree and bounded number of neighbouring vertices. 
Complexity and approximation results for the connected vertex 
cover problem on graphs and hypergraphs [30], vertex cover 
on dense hypergraphs [27], set covering on hypergraphs [28] 
and vertex cover on k-partite k-uniform hypergraphs [29] have 
been obtained. Independent set problems are closely related to 
covering problems. A comprehensive study of several approx- 
imation algorithms on hypergraph independent set problems 
can be found in [31]. However, the appUcations of EAs and 
ACO to these problems are still unexplored. 

In the present work, theoretical analysis of ACO is car- 
ried out for the minimum weight edge cover problem on 
generalized hypergraphs. The minimum vertex cover and 
the maximum weak-independent set problems are shown as 
special cases of the minimum weight edge cover problem. 
An incremental construction procedure of ACO is employed 
in this paper. The expected optimization time is calculated 
for the MM AS* algorithm with the impact of pheromone 
values and heuristic information investigated in detail. The 
expected number of iterations for finding a minimum weight 
edge cover is found to be of exponential order for generaUzed 
hypergraphs without heuristic information present. However, 
for certain families of hypergraphs, we show that with proper 
adjustment of the heuristic parameter, expected optimization 
time of constant order can be achieved. 

The rest of the paper is organized as follows. Section II 
introduces the fundamental concepts of a hypergraph and 
gives the problem definitions. Section III gives the detailed 
description of the MMAS* algorithm and the construction 
procedure for the minimum weight edge cover problem. In 
section IV, the upper bound on the expected optimization 
time is analysed. In section V, the solution methods for the 
minimum vertex cover and the maximum weak-independent 
Set problems on hypergraphs are presented as a special case 
of the minimum weight edge Cover problem. Section VI 
concludes the paper with a brief discussion on the obtained 
results. 

11. Preliminaries 

This section gives the definitions to different covering and 
independent set problems in case of hypergraphs. 

A. Basic Definitions: 

A Hypergraph H is a pair H={V, E), where 
V={v\,V2, ■■■,Vn} is a set of discrete elements known 
as vertices and £^={ei, 62, ....€„} is a collection of non- 
empty subsets of V such that IJI^i ^'i - known as 
hyperedges. Thus, a hyperedge typically contains any number 



of vertices. In the rest of the paper, the terms edge and 
hyperedge will be used interchangeably. 

The size, or the cardinality, \e\ of a hyperedge is the number 
of vertices in e. A hypergraph is known as k-uniform if all the 
hyperedges have cardinality k. A pendant vertex is a vertex 
which is contained in only one hyperedge. 

A hypergraph H(V,E) is vertex-weighted if every vertex 
in V is assigned a weight. Similarly, a hypergraph is edge- 
weighted if every edge in E is assigned a weight. In this 
paper, a weighted hypergraph will refer to an edge-weighted 
hypergraph, unless stated otherwise. 

The dual of a hypergraph, H(V,E) is a hypergraph//*, whose 
vertices and edges are interchanged with H. Thus, for each 
vertex v in H, there is an hyperedge e* in //*, formed by 
vertices in H* which correspond to hyperedges incident on v 
in H. Thus e* = {e^ G E\v € ej}. 

B. Covering Problems: 

For a hypergraph H(V,E), a Vertex Cover or a Hitting Set is 
a subset of V that intersects every edge of // in at least one 
vertex. 

Similarly an Edge Cover of a hypergraph H is a subset of 
E that contains all vertices in V. Edge Cover is a special case 
of the Set Cover problem. From the definitions, a hitting set 
in // is a set cover in the dual graph H*. Mathematically, 
Vertex Cover: A set 5 C V such that V e € E:\ef]S\ > 1. 
Edge Cover: A set J C £ such thatV v e V 3 e e T: v e e. 

The Minimum Vertex Cover and the Minimum Edge Cover 
of a hypergraph are optimization problems to find the vertex 
cover and the edge cover of minimum cardinality in a given 
hypergraph H, respectively. Similarly, the Minimum Weight 
Vertex Cover and Minimum Weight Edge Cover problems for 
edge-weighted and vertex-weighted hypergraphs respectively 
would mean to find the edge cover and the hitting set of 
minimum weight. 

C. Independent Set Problems: 

The Independent Set problem has two variants in case 
of hypergraphs. A Weak-Independent Set in // is a subset 
of V that doesn't contain any edge of H completely. If an 
independent set intersects any hyperedge in E in at most one 
element, then it is a Strong-Independent Set. Thus for graphs, 
both the above categories of independent sets are identical. In 
this paper, independent set refers only to Weak-Independent 
Set unless stated otherwise. Mathematically, 
Weak-Independent Set: A set 1^ <Z V such that V e € 
E:\e{^I^\ < \e\. 

Strong-Independent Set: A set Is <Z V such that V e € 
E:\e[^h\ < 1. 

The goal of the Maximum Weak-Independent Set of a 
given hypergraph H is to find a weak-independent set in 
H of maximum cardinality. Similarly for a vertex-weighted 
hypergraph H, the Maximum Weight Weak-Independent Set 
problem is to find a weak-independent set in H of maximum 
weight. 



III. Algorithm 

In this section, we propose an MMAS* (Max-Min Ant Sys- 
tem) ACO algorithm in the context of Minimum Weight Edge 
Cover problem in an edge-weighted hypergraph. MMAS* 
is an iterative algorithm, which in each step, constructs a 
new candidate solution and keeps track of the best-so-far 
solution. The construction procedure uses a construction graph 
in which the vertices are the elements to be chosen to 
generate a new solution. As the name suggests, the pheromone 
values are bounded between a maximum and minimum value. 
Thus a typical iteration of MMAS* consists of three steps : 
a)Constmction of a new solution; b) Selection of the generated 
solution; c) Pheromone update. The heuristic information is 
kept unchanged throughout. 



A. Construction of a new solution: 

In each iteration, a new solution is obtained by the procedure 
Construct, as described in Algorithm 1. This is an iterative 
procedure, in which the artificial ant selects the components 
of the candidate solution, one at a time, by random walk 
on the construction graph. For the edge cover problem, the 
components of the solution are the edges of the input hyper- 
graph H(V,E). Thus, we choose the construction graph C(H) 
to be a directed graph, with vertices {eo, ei, 62, ....em}, where 
CiVi > corresponds to an edge of H and eo is a starting 
vertex. The edge set U of C(H) is given as: 
U={iei, ej) I < i < m,l < j < m,i ^ j} 
i.e., in the graph C(H), there is no edge pointing to the starting 
vertex, cq. 

For the edge cover problem, we first select all the edges 
which contain at least one pendant vertex to be part of the 
solution by default as per the definition. Then, we make the 
construction graph omitting all the vertices that correspond 
to the hyperedges containing pendant vertices in the original 
hypergraph. In the first step, the artificial ant is at the starting 
node eo. In each subsequent step, it goes probabilistically to 
another node in the construction graph avoiding multiple visits 
to a particular node and only if the new node contains at 
least one new vertex of the hypergraph which is not covered 
by the previously visited nodes (i.e. hyperedges of the given 
hypergraph). Let us assume, without loss of generahty, that 
the i*'* node traversed is denoted by ej. Hence, we define the 
notion of feasible neighbourhood Nifik) for a given node 
with the previously visited nodes being {ei, 62, e/c_i, Cfe} 
as follows: 

N{ek) := E' \ {e' £ E' \ \/ V € e',v G ei U e2 U ... U efc}, 
where E' = {E\ {ei, 62, efc}). 

Thus, at every step, a new node is selected by the ant 
within its feasible neighbourhood. The selection process 
is determined by the current pheromone values r of the 
edges and the heuristic information r] of the vertices of the 
construction graph. The probability of selecting any e^+i in 
the (A; + 1)*'' step is given as: 



e'eiV(efe) 

10, 



otherwise 



where a and (3 are two positive parameters which control the 
relative strengths of pheromone value and heuristic informa- 
tion respectively in choosing the new edge of the hypergraph 
or the node on the construction graph. The construction 
procedure continues till we get a set of edges such that all 
the vertices are covered at least once. So, in each step after 
choosing an edge, we keep track of all the vertices covered so 
far. 

Since, we are interested in the minimum weight edge cover, 
we would like to prefer edges with less weight and containing 
more number of vertices. Thus, for a hyperedge e, we define 
the heuristic information as rie = -r^. 



Algorithm 1: The Procedure Construct on (C{H),T,r]) 

begin 

V = (t) 
E' = 4) 
e 4- eo 

while V Ao 

Choose a node efe+i G N{ek) after choosing Cfc 

with probabiUty 1^'°'°+^'"'°''+^ ^ 

e'eAf(efc) 

E' ^E'\J Ck+l 

V = V'yj v', W G efe+1 
_ retum E' 



B. Fitness Evaluation and Selection: 

The fitness of a candidate solution x, /(x)={ei, e2, •••efc} is 
the total weight of the edges in x. Our objective is to minimize 
this fitness value. Formally, the fitness function is stated as 
follows: f{x) := w{ei). 

The new solution replaces the best-so-far solution if its 
fitness value is strictly less than the fitness value of the best- 
so-far solution. 

C. Pheromone update: 

Initialization of pheromone values is done at the start of the 
MMAS* algorithm in such a way that their sum is equal to 
l.Thus, T(„,^) =l/|f/|; y{u,v)€U. 

After an initial solution x* is obtained using Construct, the 
pheromone values are updated using the procedure Update 
[18]. Update(T, X*) = r' is defined as follows: 



V(e,e')Gf/:r(',^,,) 



' h, if (e, e') G x* 
I, otherwise 



where h and / are two pheromone values and h > I. In 
the subsequent iterations, the procedure Update is used to 
update the pheromone values, only if the new solution obtained 
from the procedure Construct replaces the best-so-far solution 
obtained till the previous iteration. Thus, all the nodes of the 
construction graph selected in the best-so-far solution have 
pheromone value h on the incoming edges to them. The 
remaining edges have pheromone value I. 

See Algorithm 2 for the complete algorithm MMAS*. 

Algorithm 2: The Algorithm MMAS* on H(V,E) 
begin 

Tu,v^l/\U\,y{u,v)GU 

X* ^ Construct{C{H),T,r]) 
T ^ Update(T,x*) 
while true do 

X 'h- Construct{C{H),T, rj) 
if /fxj < /fx*) then 



L X* 



Update{T, X*) 



{k-i)l 



(k - i)l 



{m-k-i)h + {k- i)l ((m - k)h -\- kl) - i{h + I) 

k — i 

((m - k)cn + k)- i{c„ + 1) 



The probability that the MMAS* algorithm gives the opti- 
mal solution is at least: 



fc-i 

n 



k — i 



fc!((m — fc)c„)! 
((m - fc)c„ + fc) - i{cn + 1) ^ {{m - k)cn + ky. 



Hence, the expected optimization time is found to be of 



O 



((m — k)cn + ky. 
k\{{m-k)cny 



Now, we move to the second scenario, where we analyse 
the upper bound for the expected optimization time when 
a = and /3 = 1. Here we examine the behaviour of the 
algorithm in the presence of heuristic information. 



IV. Upper Bounds 

In this section, we show the upper bounds on the expected 
optimization time of the proposed MMAS* algorithm on the 
minimum weight edge cover problem of hypergraphs. The 
results are expressed in terms of the parameters a, P, the 
pheromone bounds h and I, the number of hyperedges m, and 
the number of hyperedges in the optimal solution k and some 
other parameters which will be introduced later in this section. 
First, we do the analysis for generalized hypergraphs and then 
we show constant order expected optimization time for certain 
famiUes of hypergraphs. 

A. Upper Bounds for the Generalized Hypergraph 

In case of generaUzed hypergraphs, we separately 
investigate the influence of pheromone values and heuristic 
information respectively. 

Theorem 1: Let a = 1 and j3 = and < I < h, and let 

h 



. Then the expected optimization time of the MMAS* 



IS 



^ ( ((m-fc)c„+fc)! \ 
^ \ ((m-fe)c„)!fe! )■ 



Proof: For the purpose of calculating the upper bound, 
we first need to find a lower bound on the probability that the 
MMAS* algorithm gives the optimal solution. We consider 
the worst possible case . Thus, we assume that there are no 
hyperedges having vertices with degree 1, the k hyperedges 
belonging to the optimal solution have pheromone level / 
and the rest of the hyperedges have pheromone level h. The 
probability of selecting a hyperedge belonging to the optimal 
solution in step i is at least: 



Theorem 2: 



Choosing a = and (3 = 1, the 



expected optimization time of the algorithm is bounded 

k 



(m-k) 



Proof: Let ei,e2, be the hyperedges of the hy- 
pergraph and r]2. -Vm be their corresponding heuristic 

information values. Let T = {r]i,r]2, , %, , Ty^jdenote 

the set that contains the heuristic information of all the hy- 
peredges. Without loss of generality, let S = {rji , 772, iVk} 

be the set that contains the heuristic information of the edges 
present in the optimal solution. Let Pi denote the probability 
that the i*'^ hyperedge selected during the course of one 
iteration, belongs to S. Also, let Ri be the set containing the 
elements of S which have not been included in the solution 
till the {i — l)th iteration. 



E 

%e((T\s)ufli 



> 



/ J Irmn 



> 



%e(T\s) 



/ ^ Irmn 



/ y 'Irmn ' / j 'Ij / j Irmn ' / j ' Imax 

Vj&Ri rije{T\S) rij&R, r/je(T\S) 



> 



/3 



1 



^min + {rn- k)r]\ 



'max 



where rjmin = min {rj \ -q e T} , rjmax = max {»7 | ?7 G T} 
and Cj = r]j — rjmin for r]j & S , j = 1 Xo k, are 



positive constants. Therefore, Pmin, the lower bound on the 
probability of choosing S as the solution is: 



1 



MIN 



\ Vmin J 



Substituting /? = 1 and taking the inverse of Pmin, we find 
the expected optimization time to be of 



O 



1 + 



Vmin 



1 k 



(m — k) 



i.e. O 



■fir, 



(to — k) 



It can thus be noted that the expected optimization time 
of the algorithm depends on the nature of the graph. In the 
worst case though, the parameter k stated above, has an upper 
bound of n. It can be shown that there exists a certain family of 
graphs for which the expected optimization time is of constant 
order. 

B. Upper Bounds for Specific Instances of Hypergraphs 

We consider the case where a = and /? = 1 and obtain 
results for hypergraphs satisfying specific properties. 

Definitions: Before proceeding to theorem 3, let us define 
the terms ijimax and rj/min- 

Vimax = max {r]\r] eT\S} 

V'min = min {r]\r] eS} 

where S and T are as defined in the proof of theorem 2. 

Theorem 3: If "qimin > Vimax (^nd we choose 



/3 > 



log{k{m — k)) 



the expected optimization time 



loairilmin/Vlma^) 

for the algorithm is of constant order 



Proof: Using the definitions mentioned in the proof of 
Theorem 2 and proceeding in a similar fashion, we obtain: 



Therefore, P'min, the lower bound on the probability of 
choosing S as the solution is: 



P' 



^+{^^)\rn-k) 



In the expression of P'min^ if we choose 
^ ^imox j 1^^^ _ bounded above by ^, we 



would obtain: 



P'min > ( 1 + 



1 

> - 

e 



Then, the expected optimization time of the algorithm 
would be bounded above by the inverse of P'^in i-®- 
0(e), which is of constant order. 

For the above assertion to be true, the condition 
< T must hold for some non- 



negative j3. Two cases may arise. 

This case will never hold good for any non-negative /3, and 
hence is rejected. 



Case 2: 



( Vlmax ^ 



P log 



< 1 



d> 



(to — k) 
< log 



< 



k 
1 



fc(m — k) 
log{k{rn — k)) 



log ( Hl^) 



Hence, it can be concluded that, if T^^mm ^ 'Himax 

and we 



E E(^'i^.n + Cj)+ E ^1 foose/3> 



log{k{m—k)) 



the expected optimization time 



> 



r,^e((T\S)Ufl« 



Vj^Ri 



> 



r)je{T\S) 



/ J I mm 

VjeRi 



/ J I mm ' / J ij / J I mm ' / j Umax 

VjeRi Vje{T\S) rijeRi r,je{T\S) 



I mm 

~ ^'Lin + k)4max 1 + {^^^^ \m - k) 

where Cj- = r]j — rj'^in for 77^ G 5 , j = 1 to A; are positive 
constants. 



log{ril_mi^/riimax) 

for the algorithm is of constant order. 



Some instances of hypergraphs for which the above 
conditions hold good and thus the algorithm has constant 
order expected optimization time, are discussed below. 

Instance 1: We shall construct a weighted complete 
r-uniform hypergraph H = {V,E) on the vertex 

set V = {1,2, ,n} and the hyperedge set 

E = {ei,e2, ,em}, where each hyperedge e is treated as 

a set of vertices. Let 5 = and T = </> be empty sets and 
i be initialized to 0. Let weight of an edge e be denoted by 
w{e). We construct the hypergraph as follows: 
(i) increment i and select a new hyperedge e^ € {E\ S) such 



that T n e ■ = (/) 

{ii) assign ^(eQ = 1. Update S = SU {e^} and T = T U 

(iii) ii i < repeat steps (i), {ii) and (iii). Else, goto 
step {iv) 

(iv) increment i and select a new hyperedge such that 

(F\r)Ce^ 

(v) assign w(e^) = 1. Update S = SU {e^} and T = T U 

set w{ej) = randj V e € {E \ S), where randj is a 
random integer > 1 

Clearly, for the above hypergraph H, the minimum weight 
edge cover is the set S. Also, since H satisfies the conditions 
of Theorem 3, expected optimization time would be of 
constant order. 

Instance 2: We construct an unweighted hypergraph 

H = {V,E) on the vertex set V = {1,2, ,n} and the 

hyperedge set E = {ei, 62, , em}- Let S = (j) and T = (j) 

be empty sets and i be initiaUzed to 0. Let p be an integer 
such that 2 < p < n. 

(i) increment i and construct a new hyperedge such that 
e- ^ 5' as follows: 

ifp< \V\T\, then C{V\T) 
else {V\T)C e\ 

{ii) update p = q such that 2 < q < p, S = S (J {ci} and 
T = TUei 

{Hi) if |r| < \V\, then repeat steps {i) and {ii), else goto 
step {iv). 

{iv) terminate construction procedure if so desired or if 
the maximum limit for the number of hyperedges i.e. 2" 
is reached, {v) update p = q such that 2 < q < p and 

5 = 5U{eJ . 

{vi) increment i and construct a new hyperedge e- such that 

e'^^ S and \ei\ = p. 

The expected optimization time of the algorithm for the 
above mentioned Instance 2 of hypergraphs, is of constant 
order. 

Proof: Let the set of all hyperedges belonging to S 
that were constructed using step (?) be denoted as ^. Let 
$ = S \ 'i'. Let the last hyperedge that was constructed 
using step (i) be denoted by eiast- Since the hypergraph is 
unweighted, we may assume that each hyperedge has weight 
equal to 1. Now, in order to show that is the minimum 
edge cover of the hypergraph H, we need to show that no 
subset of ^ can be replaced by any subset of $ while not 
increasing the total number of edges used in the edge cover 
Let U C {^ \ {eiast}) and V C ^. For V to be able to 
replace U, there must exist some V s.t. \V\ < \U\ and V 
covers all vertices covered by U. However if \V\ < \U\, the 
number of vertices covered by V is less than the number 
of vertices covered by U since V(e[/,ey) s.t. eu & U and 
ev £ V, \eu\ > |ey| and elements of U are disjoint while 
elements of V are not necessarily disjoint. 



Hence, for no U can a corresponding V exist, which can 
replace U while satisfying the condition |V"| < |C/|. (1) 

Also, for no {U U {eiast}) can a corresponding V exist, 
which can replace {UL){eiast}) while satisfying the condition 
\V\ < \U\. (2) 

Hence, using statements (1) and (2), it can be concluded 
that, no subset of can be replaced by any subset of $ while 
not increasing the total number of edges used in the edge 
cover. Thus, Vl/ is a minimum vertex cover. 

Applying theorem 3, there exists 13 s.t. the expected opti- 
mization time of the algorithm is of constant order. 

■ 

V. Other Hypergraph covering problems 

A. Minimum Vertex Cover Problem 

The vertex cover, as defined earlier, is a set of vertices such 
that each hyperedge contains at least one vertex belonging to 
this set. The minimum vertex cover problem is finding a vertex 
cover of minimum cardinality. 

The minimum edge cover problem is a special case of the 
set cover problem where we treat each set as a hyperedge and 
the union of all the given sets as the set of vertices. Similarly, 
the minimum vertex cover problem is a special case of the 
hitting set problem. It is known that the set cover problem 
and the hitting set problem are equivalent reformulations of 
one another. Thus, for hypergraphs, the minimum vertex cover 
problem is an equivalent reformulation of the minimum edge 
cover problem. 

Hence, we can address the minimum vertex cover problem 
by reducing to the minimum edge cover problem, for which 
we have already provided an MM AS* algorithm. For this 
purpose, we take the given hypergraph H{V, E) as input, find 
its dual hypergraph H* {V* , E* ) with unweighted edges and 
then apply the MMAS* algorithm for finding the minimum 
edge cover to H*. The MMAS* algorithm gives a set of 
hyperedges as output which is a subset of E* . By the definition 
of dual of a hypergraph, these set of hyperedges of H* given 
as output, are vertices of the hypergraph H and hence form 
a solution for the minimum vertex cover problem for the 
hypergraph H. 

B. Maximum-Weak Independent Set Problem 

We know that a given set of vertices forms a weak- 
independent set if and only if its complement is a vertex cover. 
Hence, the solution to the weak independent set problem can 
be obtained directly from the solution of the minimum vertex 
cover problem, by taking the complement of the minimum 
vertex cover 

It should also be noted that the maximum chque problem 
can solved by taking a complement hypergraph H' of the 
input hypergraph H and then finding the maximum weak- 
independent set for H' . 



VI. Conclusion 

In the present work, we obtain results for the upper bound 
of expected optimization time of MM AS* ACO algorithm 
for the minimum weight edge cover problem on generalized 
hypergraphs by taking into account pheromone values and 
heuristic information separately. The results obtained are of 
exponential order with respect to the number of hyperedges 
for generalized hypergraphs. However, for certain families 
of hypergraphs, by suitably choosing the heuristic parameter, 
expected optimization time of constant order could also be 
obtained. We show two such theoretical instances where the 
above holds true. Also, we show that the minimum vertex 
cover and maximum weak-independent set problems are spe- 
cial cases of the minimum weight edge cover problem. 

Further work may be taken up in the following directions. 

• Investigation of different ACO construction procedures 
and corresponding runtime analysis for the covering prob- 
lems on specific classes of hypergraphs, i.e. k-uniform, 
k-partite, dense/sparse hypergraphs, etc. 

• In-depth analysis of the impact of pheromone values and 
the pheromone evaporation factor on the running time. 

• Analysis of lower bounds on the expected optimization 
time for different classes of hypergraphs. 

• Experimental investigation of the algorithm's perfor- 
mance. 

• The combined effect of pheromone and heuristic infor- 
mation on performance of the algorithm. 

• Extension of the present work to other NP-Hard/NP- 
Complete problems on hypergraphs, eg. coloring prob- 
lems, minimum spanning tree problem, etc. 

The present study thus points to several interesting open 
areas in the theoretical ACO research that may contribute to 
a better understanding of the runtime behaviour of ACO for 
computationally harder problems. 

References 

[1] M. Dorigo, V. Maniezzo, and A. Colorni, "The ant system: An autocat- 
alytic optimizing process", Dipartimento di Elettronica, Politecnico di 
Milano, Milan, Italy, Tech. Rep. 91-016, 1991. 

[2] M. Dorigo and T. Stiitzle. Ant Colony Optimization. MIT Press, Cam- 
brigde, MA, 2004. 

[3] B. DoeiT, F. Neumann, D. Sudholt, and C. Witt. "On the runtime analysis 

of the 1-ANT ACO algorithm". In Proc. Genetic and evolutionary 

computation (GECCO '07j,London. ACM,pp. 33-40,2007. 
[4] T. Stiitzle and H. H. Hoos, "MAX-MIN ant system" , Future Gener. 

Comput. Syst., vol. 16, no. 8, pp.889-914, 2000. 
[5] T. Stiitzle and M. Dorigo, "A short convergence proof for a class of ACO 

algorithms", IEEE Transactions on Evolutionary Computation, vol. 6, no. 

4, pp. 358-365,2002. 
[6] W. J. Gutjahr,"A graph-based ant system and its convergence". Future 

Gener. Comput. Syst., vol. 16, no. 9, pp. 873-888, 2000. 
[7] S. Droste, T. Jansen, and I. Wegener. "On the analysis of the (1+1) 

evolutionary algorithm". Theoretical Computer Science, vol. 276, no. 

l-2,pp.51-8I, 2002. 
[8] M. Dorigo and C. Blum, "Ant colony optimization theory: A survey". 

Theoretical Computer Science, vol. 344, no. 2—3, pp. 243—278, 2005. 
[9] W.J. Gutjahr, "Mathematical runtime analysis of ACO algorithms: Survey 

on an emerging issue". Swarm Intelligence vol. l,no. 1, pp. 59— 79 (2007). 



[10] B. Doerr,D.Johannsen,"Refined runtime analysis of a basic ant colony 
optimization algorithm". In Proc. IEEE Congr Evol. Comput., Piscat- 
away, NJ: IEEE Press,pp. 501-507,2007. 

[II] F. Neumann and C. Witt, "Runtime analysis of a simple ant colony 
optimization algorithm". In Proc. 1 7th Int. Symp. Algorithms Com- 
p«r.(/SAACj,Kolkata, India. LNCS vol. 4288, pp.618-627,Beriin, Ger- 
many: Springer- Verlag, 2006. 

[12] W. J. Gutjahr, "First steps to the runtime complexity analysis of ant 
colony optimization", CompMf. Oper. Res., vol. 35, no. 9, pp. nil — 
2727,2008. 

[13] F. Neumann, D. Sudholt, and C. Witt,"Comparing variants of MMAS 
ACO algo- rithms on pseudo-boolean functions". In Proc. Stoch. Local 
Search Algorithms (SLS'07), LNCS 4638. Brussels, Belgium: Springer- 
Verlag, pp. 61-75,2007. 

[14] W. J. Gutjahr and G. Sebastiani,"Runtime analysis of ant colony opti- 
mization with best-so-far reinforcement" Methodology and Computing in 
Applied Probability, vol. 10,no 3,pp.409-433, 2008. 

[15] F. Neumann and C. Witt, "Ant colony optimization and the minimum 
spanning tree problem", in Proc. Learn. Intell. Optimization-LION 2, 
LNCS vol. 5313, Berlin, Germany: Springer- Verlag, pp. 153-166,2007. 

[16] N. Attiratanasunthron and J. Fakcharoenphol, "A mnning time analysis 
of an ant colony optimization algorithm for shortest paths in directed 
acyclic graphs", /n/orm. Process. Lett., vol. 105, no. 3, pp. 88—92, 2008. 

[17] C. Horoba and D. Sudholt. "Running time analysis of ACO systems 
for shortest path problems". In Proc. Stoch. Local Search Algorithms 
(SLS'09),LNCS vol. 5752, pp. 76-91. Springer,Heidelberg, 2009. 

[18] T.Kotzing, P.K. Lehre, F. Neumann, and PS. 01iveto,"Ant colony 
optimization and the minimum cut problem". In Proc. Genetic 
and evolutionary computation (GECCO 'iO),Portland,USA. ACM: 
pp. 1393- 1400,2010. 

[19] Y. Zhou. "Runtime analysis of an ant colony optimization algorithm 
for TSP instances", IEEE Transactions on Evolutionary Computation, yo\. 
13,no. 5,pp. 1083-1092, 2009. 

[20] T. Kotzing, F. Neumann, H. Rglin, and C. Witt, "Theoretical properties 
of two ACO approaches for the traveling salesman problem". In Proc. of 
ANTS'2010,UiCS vol. 6234, pp. 324-335, 2010. 

[21] D. Sudholt. Computational Complexity of Evolutionary Algorithms, 
Hybridizations, and Swarm Intelligence. PhD thesis, Technische Universit 
at Dortmund, 2008. 

[22] F. Neumann and C. Witt.Bioinspired Computation in Combinatorial 
0/?H'nifea?io«. Springer- verlag,Berlin,Germany.2010. 

[23] R.M. Karp, "Reducibility among combinatorial problems". Complexity 
of Computer Computations, Plenum Press,pp. 85 — 103, 1972. 

[24] N. Bansal and S. Khot, "Inapproximability of hypergraph vertex cover 
and applications to scheduling problems". In Proc. of ICALP' 1 0,lJ^CS 
6 198,Springer- Verlag, Heidelberg, pp. 250-261,2010. 

[25] M. Okun. "On approximation of the vertex cover problem in hyper- 
graphs". Discrete Optimization, vol.2,no. 1, pp.101 — 111, 2005. 

[26] I. Dinur, V. Gurus wami, S. Khot, and O. Regev. "A new multilayered 
PCP and the hardness of hypergraph vertex cover", /« Proc. 35th ACM 
STOC, pages 595-601, 2003. 

[27] J. Cardinal, M.Karpinski, R.Schmied, C.Viehmann, "Approximating Ver- 
tex Cover in Dense Hypergraphs", arXiv:10I2.2573vl [cs.DS].2010. 

[28] M. Krivelevich, "Approximate set covering in uniform hypergraphs", / 
Algorithms, m\. 25,no. I,pp. 118-143,1997. 

[29] V. Guruswami and R. Saket,"On the Inapproximability of Ver- 
tex Cover on k-Partite k-Uniform Hypergraphs", /«. proc. ICALP'IO 
LNCS,Springer-verlag,Heidelberg, vol. 6198 , pp. 360-371,2010. 

[30] B. Escoffier.L. Gourves.J. Monnot, "Complexity and approximation re- 
sults for the connected vertex cover problem in graphs and hyper- 
graphs" , Journal of Discrete Algorithms, vol 8, no l,pp 36—49,2010. 

[31] E. Losievskaja,Approximci;;o« Algorithms for Independent Set Problems 
on IIypergraphs,Vhd thesis. The School of Computer Science, Reykjavk 
University.2009. 



